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ABSTRACT 

Current cosmological data indicate that our universe contains a substantial component of dark vacuum 
energy that is driving the cosmos to accelerate. We examine the immediate and longer term consequences 
of this dark energy (assumed here to have a constant density) . Using analytic calculations and supporting 
numerical simulations, we present criteria for test bodies to remain bound to existing structures. We 
show that collapsed halos become spatially isolated and dynamically relax to a particular density profile 
with logarithmic slope steeper than —3 at radii beyond T2oo- The asymptotic form of the space-time 
metric is then specified. We develop this scenario further by determining the effects of the accelerating 
expansion on the background radiation fields and individual particles. In an appendix, we generalize 
these results to include quintessence. 

Subject headings: cosmology:theory — large-scale structure of universe — dark matter — dark energy 
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1. INTRODUCTION 

The past several years have witnessed an impressive so- 
lidification of our estimates for the basic cosmological pa- 
rameters. Measurements of the cosmic microwave back- 
ground radiation indicate that the cosmos is spatially flat 
(e.g., Hanany et al. 2000). Observations of Type la super- 
novae suggest that the universe is accelerating (Riess et 
al. 1998; Perlmutter et al. 1998; Garnavich et al. 1998). 
These data, in conjunction with additional observations 
of the large-scale galaxy distribution (e.g., Peacock et al. 
2001) and the Hubble constant (Freedman et al. 2001), ar- 
gue for a specific cosmological model with matter density 
il mj o — 0.3, dark vacuum energy density fl v .o ~ 0.7, cur- 
vature constant k = 0, and Hubble constant Ho —70 km 
s" 1 Mpc -1 . This cosmological model has recently received 
impressive validation by results from the WMAP satellite 
(Bennett et al. 2003). In combination with data on the 
Hubble constant, large-scale structure, and they Lyman-a 
forest, WMAP constrains the values of these parameters 
to a few percent (Spergel et al. 2003). 

In this newly consolidated cosmological model, the most 
unexpected property is the large (apparent) energy den- 
sity in the form of dark vacuum energy. The substantial 
value of Q v fi = 0.7 produces a correspondingly large effect 
on past cosmological evolution. In particular, the age of 
the universe is somewhat longer than in a flat cosmology 
with no vacuum energy and more consistent with other age 
indicators (e.g., Carroll, Press, & Turner 1992). But, by 
far, the most striking consequence of this vacuum energy 
lies in our cosmological future. 

If the universe is already starting to accelerate, as indi- 
cated by the observationally implied value Q V) o = 0.7, then 
structure formation is virtually finished. In the relatively 
near future, the universe will approach a state of expo- 
nential expansion and growing cosmological perturbations 
will freeze out on all scales. Existing structures will grow 
isolated. In the face of such desolation, we would like to 



know more quantitatively the conditions required for the 
formation (collapse) of future cosmological structures and 
the conditions required for small bodies to remain bound 
to existing structures. We would also like to know the 
asymptotic form of the existing cosmological structures, 
in particular the dark matter halos. By answering these 
questions, we can determine how much structure forma- 
tion remains to occur in the future of our universe. 

Given the relatively well constrained parameters of our 
universe, the future evolution of cosmological structure can 
now be predicted with some confidence. Indeed, several 
recent papers have begun to explore this issue. Possible 
future effects of vacuum energy density were outlined in 
a recent review of our cosmic future (Adams & Laughlin 
1997, hereafter AL97). As the universe accelerates, cur- 
rently visible galaxies are redshifted out of view and will 
become inaccessible to future astronomers (Loeb 2002). 
Simulations of future structure formation have been done 
for the case of a cosmological constant with a focus on our 
local portion of the universe (Nagamine & Loeb 2002). 
Similar issues have been explored semi-analytically, in- 
cluding basic effects of quintessence (Chieuh & He 2002; 
Gudmundsson & Bjornsson 2002). A comprehensive list 
of papers related to the future evolution of the universe is 
compiled in Cirkovic (2003). 

In this paper, we consider the future evolution of struc- 
ture formation with a constant density of dark vacuum 
energy. We extend previous work by deriving analytic 
estimates for the conditions required for the collapse of 
structures, and by analyzing the results from a suite of 
numerical simulations of future structure formation. The 
future evolution of cosmic structure can be viewed in two 
related ways. The traditional approach considers whether 
or not a given region of the universe with overdensity 
So = (p— (p))/(p) will collapse. Given the strong suppres- 
sion of future structure formation, however, relatively little 
will happen in terms of new formation. A related question 
is to ask whether or not test bodies will remain bound to 
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existing structures. This issue operates on a wide range 
of spatial and mass scales: Will the local group remain 
bound to the Virgo cluster? Will a satellite dwarf galaxy 
remain bound to the Milky Way? What bodies would 
remain bound to an isolated star in the face of acceler- 
ated expansion? This paper develops both approaches — 
analytically in §2 and numerically in §3 — and elucidates 
the relationship between them. Our numerical simulations 
also show that the density profiles of dark matter halos ap- 
proach a nearly universal form. Every dark matter halo 
grows asymptotically isolated and becomes the center of 
its own island universe; each of these regions of space-time 
then approaches a universal geometry. In §4, we determine 
the sphere of influence of existing structures and find the 
future time at which they grow cosmologically isolated. 
Next, we consider the implications of cosmic acceleration 
for the background radiation fields, suppression of particle 
annihilation, and the long term geometry of space-time. 
Because the equation of state of the dark vacuum energy 
is not completely determined, we generalize these results 
(in an Appendix) to consider dark energy that varies with 
time (equivalently, the scale factor). 

2. ANALYTICAL DESCRIPTIONS OF STRUCTURE 
FORMATION 

Before performing detailed numerical simulations of fu- 
ture structure formation, it is useful to develop simple an- 
alytical estimates. Such results have been developed previ- 
ously to account for past structure formation in an accel- 
erating universe (see especially Lokas & Hoffman 2001). 
While most past work has emphasized the evolution of 
structure up to the present epoch, we focus here on its 
evolution into the future. Throughout this treatment, we 
assume a spatially flat universe, which at the present epoch 

has ri m ,o — 0-3 and ^ v ,o = 0.7 (flm,o + ^v,o = 1)- In the 
future, the values of fi m and v vary, but their sum con- 
tinues to equal unity. 

2.1. Collapse of over 'dense regions 

The evolution of a given region of the universe is de- 
scribed by an energy equation. For a spherical patch of 
physical size r, this energy equation can be written in the 
form 

^-^-^rX„ = £, (1) 

where the energy E is given by 
1 



E 



Hlrl[l - n m , (l + 5 ) - fi v ,o] 



(2) 



This set of equations implicitly assumes that r = Hqtq, 
i.e., that the particles are traveling along with the unper- 
turbed Hubble flow at the present epoch. Later in this 
section, we generalize this analysis to consider the case in 
which particles have already slowed down relative to the 
Hubble flow due to the past action of gravity. 

If we define the dimcnsionless variable £ = r/r and the 
dimensionless time t = Hot, then the energy equation can 
be rewritten in the simpler form 

1 + n m , (| -!)(! + So) + «v,o(£ 2 - 1) ■ (3) 

If a cosmological structure is slated to collapse in the fu- 
ture, then the effective velocity (the time derivative of £) 



K dT J 



must change sign, which requires the right hand side of 
the above equation to vanish. This requirement results in 
a cubic equation of the form 

e 3 +[^-l-§^(l + 5o)]? + §^(l + M=0. (4) 

"v,0 s 'v,0 " "v,0 

For a given (flat) cosmology (a given value of Q m fl, which 
in turn specifies f2 v .o — 1 — ^m,o), and a given overdensity 
So, the cubic has three real roots if the following constraint 
is satisfied: 

[§^(l + ^o)] 2 <^[l + ^(l + ^o)-^] 3 . (5) 

"v,0 ^ I "v,0 "v,0 

The minimum overdensity So occurs when the above con- 
straint is saturated (i.e., at equality). For a given value 
of O mi o, the vacuum energy density O v ,o = 1 — ^m,o is 
specified, and the equation has a given root. For the cur- 
rently favored cosmological model with Sl m o = 0.3, the 
root occurs for So = 17.6. 

2.2. Criterion for being bound to existing structures 

A related issue is to ask whether a small mass or test 
particle will be bound to currently existing cosmological 
structures. To carry out this calculation, we consider an 
existing object of mass M Q bj, which could be (the dark 
matter halo encompassing) a galaxy or a cluster of galax- 
ies. We then ask whether test particles — much smaller 
structures exterior to the system — will be bound to M bj 
or not. The test particles start with a radial distance ro 
at the present epoch. For simplicity, we assume that the 
distance r is much larger than the size of the collapsed 
object so that the potential of a point mass provides a 
good approximation. 

The calculation is analogous to that of the previous sec- 
tion, with the mass of the potentially collapsing system re- 
placed by the mass M ob j of the pre-existing structure. The 
energy equation can again be written in non-dimensional 
form 

(§) 2 = 1 + (<W + 0) [| - 1] + Ov,0(£ 2 - 1) , (6) 

where we have defined 
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(7) 



The parameter (3 thus measures the effective "strength" 
of the galaxy or cluster. The requirement that equation 
(6) have a turnaround point leads, as before, to a cubic 
constraint, which now takes the form 



fVo£ 3 -/3£+(fVo + /?) = 0, 



which in turn requires 



27 

/3 3 > —(n m ,o + pfn v ,o- 



(8) 



(9) 



The minimum value of the strength parameter occurs 
when the equality is saturated and we denote the corre- 
sponding value of (3 by (3* . Solving equation (9), we find 
the value f3* 5.3. We can thus determine the condition 
that must be met in order for a test body to remain bound 
to an object of mass M ob j, i.e., 



2GM obj > (3*H 2 rl . 



(10) 
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Inserting numerical values and scaling the result, we thus 
obtain the following criterion 

where we have defined /170 = H o /(70 km s _1 Mpc -1 ). 

As an immediate application of the above result (eq. 
[11]), we can determine whether or not the Milky Way 
(Local Group) will remain bound to the Virgo cluster in 
our cosmological future. The mass of Virgo is estimated 
to be Mvirgo = 5 x 10 13 — 10 14 M Q and its current distance 
from the Milky Way is about ro = 16 Mpc (e.g., Jacoby 
et al. 1992). The requirement specified by equation (11) 
is not met — the mass falls short by a factor of 100 — so 
that the Milky Way is not destined to be bound to Virgo. 
This result is verified by numerical simulations (see §2 and 
Nagamine & Loeb 2002). 

Equation (11) thus defines an effective sphere of influ- 
ence for any given astronomical object — when considered 
as isolated in a background universe dominated by a cos- 
mological constant. For our Milky Way galaxy, this sphere 
has radius ro ~ 0.7 Mpc. For an isolated star (i.e., a free- 
floating star not associated with a galaxy), with typical 
stellar mass = 0.5 M Q , the sphere of influence has size 
ro « 55 pc. This size scale suggests that isolated binary 
star systems can remain safely bound. Although isolated 
pairs of galaxies can also remain gravitational bound, they 
live much closer to the brink of instability. 

The above analysis finds the sphere of gravitational in- 
fluence for test bodies that arc moving along with the 
Hubble flow at the present epoch. However, due to the 
past evolution of the universe, bodies that are now out- 
side galaxies (or clusters) can be slowed down relative to 
the Hubble flow due to the action of gravity in the past. 
In particular, the particles simulated in N-body simula- 
tions can be slowed down relative to the Hubble flow. To 
include this effect in the analysis, we modify equation (6) 
to take the new form 

(§) 2 = A + ("m,o +/?)(£-!)+ fVo(£ 2 - 1) . (12) 

where the constant A represents the fact that the test par- 
ticles have been slowed down relative to the Hubble flow 
(and rest of the quantities are the same as in eq. [6]). 
The value A = 1 corresponds to particles moving with the 
Hubble flow, so that particles that have been slowed will 
have A < 1. The case A = provides a benchmark, where 
test bodies have zero velocity (are turning around) at the 
present epoch. 

The requirement that equation (12) have a turnaround 
solution implies a modified form of equation (9), namely, 

27 

(/3 + 1-A) 3 > _tt v ,o(a»,o + /3) 2 . ( 13 ) 
By solving the above cubic for a given A when the in- 
equality is saturated, we find the root (3* (A) that can be 
used to define the sphere of influence of a given cluster 
or galaxy according to equation (10). For the benchmark 
case, A = 0, the root [3* « 1.1, and hence the effective 
sphere of influence is larger than the previous case by a 
factor of 1.7. 

In the following section, we compare these analytic pre- 
dictions to the results of numerical simulations. If we in- 
sert test bodies at the present epoch — at rest with respect 

2 http://www.MPA-Garching.MPG.DE/gadget/ 



to the Hubble flow — then the sphere of gravitational in- 
fluence of existing structures is described accurately by 
equation (11). In addition, although mass particles that 
begin at rest (with respect to the Hubble flow) in the dis- 
tant past are slowed down and have values A < 1 at the 
present epoch, the prediction of equation (11) works re- 
spectably well (see the following section). 



3. NUMERICAL SIMULATIONS OF STRUCTURE FORMATION 

To evaluate the analytic results of the previous section, 
particularly equation (11) defining the gravitational sphere 
for influence for an object of mass M, we have run a series 
of numerical simulations of the evolution of structure in 
a A-dominated universe into the future. All simulations 
were run using the GADGET code (Springcl, Yoshida, & 
White 2001) 2 on a parallel computer cluster at Michigan's 
Center for Advanced Computing. A flat ACDM model 
is assumed, with matter density fl m ,o = 0.3, power spec- 
trum normalization erg = 1.0, and h = 0.7, where the 
Hubble parameter Ho — 100/ikms" 1 Mpc -1 . The simu- 
lations followed the evolution of a 256.48 /i _1 Mpc region 
with 128 3 dark matter particles of mass 6.70 x 10 n h~ 1 M Q , 
starting from a = 0.258 (z = 20.8), through the present, 
and forward to a = 100. A gravitational force softening 
of 200/i _1 kpc (in physical units) was used throughout the 
computation. 

In addition to the simulation described above, which was 
evolved to a=100, we performed a similar simulation that 
evolved the a — 1 configuration forward with an added 
set of initially stationary test particles. In this run, 96 
particles were placed in a spherically symmetric fashion 
around the centers of 19 dark matter halos. In each case, 
they were placed at rest on one of 12 concentric spheres 
centered at the most bound position of the selected group. 
This simulation was also evolved to a = 100. 

The top two panels of Figure 1 show structure in a fixed 
comoving region 128 Mpc wide by 25 h' 1 Mpc deep 
at (a) the present epoch and (b) a = 100. The large-scale 
pattern of the cosmic web is well established by a = 1 
and evolves little thereafter (as emphasized by Nagamine 
& Loeb 2002). Clusters of fixed physical size shrink as 1/a 
in the comoving frame, so the halo population effectively 
condenses into a sea of "droplets" embedded in the frozen 
linear modes that define the filaments, walls, and voids 
that characterize the cosmic structure today. 

Panels (c) and (d) change perspective by showing how 
the physical region of panel (a) appears at a = 11 and 
a = 100, respectively. The physical separation between 
bound structures grows exponentially in time during the 
dcSitter expansion phase of the dark energy dominated 
era. The future is increasingly lonely. 

To compute the size of a bound halo — its gravitational 
sphere of influence defined in the previous section — we 
identify the smallest radius at which the local mean radial 
velocity is significantly larger than zero. We measure such 
sizes using two different tracers: the simulation particles 
themselves or the test particles of the second simulation. 
Since the latter are embedded at rest at a = 1, their fu- 
ture evolution should better follow, within the limits pro- 
vided by a monopole description of gravity, the prediction 
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Fig. 1. — Evolution of structure in a A-dominatcd universe. Panels a) and b) show snapshots of a comoving region of the universe 128 h^ 1 
Mpc on a side and 25 h^ 1 Mpc thick today (a = 1) and at cosmic age 91.5 Gyr in the future (a = 100), respectively. Panels c) and d) show 
regions of the same physical size as that shown in panel a) at epochs a = 11 and a = 100, respectively. The box in panel a) locates the region 
shown in b), and the box in panel c) locates the region shown in d). 
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Fig. 2. — Gravitational sphere of influence of a cosmological 
object, as a function of the mass enclosed within the sphere of in- 
fluence at a=l. The lower solid line marks the analytic prediction 
for objects starting at rest with respect to the Hubble flow (eq. 
[11]); the upper dashed line marks the analytic prediction for ob- 
jects that are at rest relative to the cluster at the present epoch 
(i.e., objects that arc marginally separated from the Hubble flow 
— see eq. [12]). These analytic predictions are compared with 
simulations, both for cluster particles (open squares) and for the 
test particles placed in the simulation at rest (stars). 



of equation (11) which assumes a pure, spherical Hubble 
flow. 

Measurements of the gravitational sphere of influence 
are shown in Figure 2. Note that the relevant mass is the 
mass within the sphere of influence at the present epoch; 
this is typically a factor of 2-3 larger than M 2 oo- The rela- 
tion between this mass and its sphere of influence for test 
particles placed in the simulation at rest at a — 0, shown 
as stars, is in excellent agreement with the predictions of 
equation (11), shown as the lower line in the figure. Val- 
ues for the simulation particles, shown as open squares, are 
slightly larger (by about 10%) because radial infall in the 
weakly non-linear regime reduces the kinetic energy asso- 
ciated with Hubble flow. These sizes are bounded from 
above by the modified relation, equation (12), with A = 0, 
shown as the upper dashed line in the figure. 

3.1. Halo Phase-space Structure 

In the far future, the character of halos' radial phase- 
space structure is markedly different from that of today. 
The upper row of Figure 3 shows the radial velocity pat- 
tern of the most massive halo at a = 1, 11 and 100. The 
halo identified at present remains the most massive at all 
future times. At each epoch, lengths are expressed in units 
of r200) the radius within which the mean enclosed density 
is 200 times the critical value. Physical velocities are ex- 
pressed in units of v 2 oo = ^/GM 2 oo A200: where M 2 oo is 
the mass within r 2 oo- The halo's physical size grows from 
M 200 = 1.08 x IO^/i^Mq at a = 1 to 2.35 x 1O 15 /i- 1 M 



a=l a=ll a=100 




Fig. 3. — Radial velocity relative to the cluster center as a func- 
tion of distance for particles in the largest halo at the epochs in- 
dicated. The upper panel shows all matter while the lower panel 
shows only those particles that lie within 2.5r200 at a = 1. 



at a = 11, and remains nearly constant thereafter. 

At a = 1, positive velocity particles at radii r/r 20 o ~ 
1 — 2 represents material that has penetrated the halo core 
and is streaming outward to large radii. This "processed" 
material mixes at these radii with a fresh stream infalling 
at v r ~ — u 20 o- The stream extends outward through a 
zero- velocity surface at r ~ 5r 20 o, and asymptotically ap- 
proaches the Hubble flow at large radii. The present-epoch 
interior phase structure displays a lumpy morphology, ev- 
idence of recent merger activity that has not yet dynami- 
cally relaxed. The infall pattern shows spikes of enhanced 
dispersion, the signatures of neighboring halos. Those 
within the gravitational sphere of influence are destined 
to merge with the central halo. 

By a = 11, phase mixing has smoothed the interior 
structure considerably, but the remains of a recent merger 
can be seen in a small tail of outflow centered at r /r 2 oo ~ 6 
(circled in the panel). The infall regime is much quieter, 
with only a few collapsed structures evident in the tail 
of Hubble flow. At a = 100, the interior phase struc- 
ture is extremely homogeneous and the infall regime is 
essentially silent; no collapsed structures perturb the flow 
within 50r 20 o- 

The lower panels of Figure 3 show how the non-linear 
material of the present cluster evolves into the future. 
Only those particles lying within 2.5r 2 oo at a = 1 are 
shown. This radius is chosen to encompass all the "pro- 
cessed" halo material at the present epoch. Although this 
boundary extends beyond most common choices for the 
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"virial radius" of a halo (Evrard, Metzler & Navarro 1996; 
Eke, Cole & Frenk 1996; White 2001), it is evident that 
some material currently within this radius is destined for 
future escape. A thin tail forms of material lifted off the 
halo, representing 0.8 and 2.6 percent (at a — 11 and 100, 
respectively) of the mass within 2.5r 20 o at a = 1. Note 
that the merger system at a = 11 circled in the upper 
panel is not as strongly evident in the lower panel, indi- 
cating that it lies beyond 2.5r 2 oo at the present epoch. 

The simple, smooth phase structure at late times sug- 
gests a long-term equilibrium. In the next subsection, we 
address the question of the eventual shape of the radial 
mass profiles of halos. 

3.2. Halo Density Profile 

Numerical studies have revealed that the non-linear den- 
sity structure of halos formed through gravitational clus- 
tering takes on a common form. The first such studies 
(Navarro, Frenk & White 1996; 1997) showed that this 
density profile can be written in the form 

4p s 



p(r) 



r/r s [l + (r/r s )f 



(14) 



where r s is a characteristic radius where the logarithmic 
slope of the density profile is isothermal: dlnp/dlnr = —2. 
The ratio r 20 o/r s defines the concentration parameter c. 
Note we have chosen the convention p s = p(r s ); in this 
case the inner density can be directly related to the con- 
centration of an NFW profile via 

1 A c 3 

Ps = T^cPc = T^r-Tj— \ TTT- sPc- (15) 

4 12 ln(l + c) — c/(l + c) 

Previous fits to this density profile have been limited to 
radii r £ r 2 oo and epochs a < 1. We present here an ex- 
tension of this form to larger radii and future epochs. 

Figure 4 shows radial density profiles derived from stack- 
ing the 50 most massive halos at each epoch displayed. 
Solid lines in the figure show binned profiles for five epochs 
from a = 1 100, while dashed lines show fits to the form 



p(r) 



Ap s 



r/r s {l + (r/r s )P}3/2 



[1 + r/r 



il+3p/2 



(16) 



Here r s is a scale radius similar to the that of equation 
(14), Too is an asymptotic radius, p is a free parameter and 
A = 2 3 p/ 2 /[1 + r s /r 00 ] 1+3 P/ 2 . We find that values p = 1.8, 



0.50 and r^ = 4.7r 2 ooa' 



6/(3p+2) 



provide fits that are 



accurate to {(Sp/ p) 2 ) 1 ^ 2 ~ 35% over the full range of radii 
and epochs examined. The scaling of r^ with expansion 
factor ensures that the profile approaches the mean mass 
density of the universe as r — > 00. 

The profile of equation (16) is steeper than the NFW 
form at radii beyond r 20 o- As r^, — > 00, the logarith- 
mic slope of the profile well beyond r 2 oo approaches 3.7 
(shown by the dot-dashed curve in Figure 4), steeper than 
the slope of 3 from the NFW case. This difference in slope 
keeps the enclosed mass from being logarithmic divergent, 
as implied by a formal extrapolation of the NFW form 
(compare eq. [14] with [16]). 

In fact, the mass of a halo in the far future will be sim- 
pler to define than it is today. At present, radial infall 
and incomplete dynamical relaxation make the choice of 
the "edge" of a cluster somewhat arbitrary (White 2001; 



Evrard & Gioia 2002). In the relatively near future, how- 
ever, halos evolve toward an equilibrium configuration that 
is bounded by an increasingly sharp zero-velocity surface 
(Figure 3). Ultimately, a meaningful and unique defini- 
tion of mass emerges, namely, all of the matter lying inte- 
rior to this well-defined zero-velocity surface. In addition, 
as shown by Figure 4, the density profile attains a well- 
defined form. 
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FlG. 4. — The asymptotic form for the density distribution of dark 
matter halos. The solid curves show the (nearly) universal form for 
dark matter halos, for a variety of epochs (from top to bottom, a = 
1, 3.38, 11.4, 38.6, and 100). Starting just after the present time, 
the dark matter halo profiles show essentially the same form, only 
the outer boundary is stretched to match onto the ever-lower den- 
sity of the background universe. The dashed curves show the fit to 
the numerical results, equation (16), evaluated at each epoch with 
the parameters given in the text. The dot-dashed curve shows the 
asymptotic form of the density profile in the limit t — > 00. 



4. LONG-TERM RAMIFICATIONS OF COSMIC 
ACCELERATION 

The considerations of the previous sections outline the 
requirements necessary for future structure formation and 
the criteria for test bodies to remain bound to existing as- 
tronomical structures. With these results in place, we can 
fill out the picture of the future evolution of our cosmos. 
In particular, we can determine the time scales for vari- 
ous bound structures to grow isolated, the corresponding 
effects on the background radiation fields in the universe, 
and the freezing out of particle annihilation processes. 

4.1. Isolation of bound structures 

Because our universe contains a dark vacuum compo- 
nent, it has a well-defined horizon scale. For unbound 
objects living in this accelerating universe, the next issue 
is thus to determine when the objects leave the horizon. 
The background universe can be described by its line ele- 
ment, which can be written in advanced time coordinates 
in the form 

ds 2 = - [1 - x V] dv 2 + 2dvdr + r 2 dfl 2 . 



(17) 
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This line element ignores the curvature due to 
gravitationally-condensed structures (see the following 
subsection). The parameter \ is related to the magnitude 
of the cosmological constant and is defined by the relation 
(in natural units) 

'27T 3 \l/2 A 2 

(18) 



x= Vl5/ 



where A is the effective temperature scale of the cosmolog- 
ical vacuum energy (A w 0.003 eV « 35 K for the presently 
suspected cosmological constant). In such a universe, the 
horizon distance r# is given by 

c 2 / 15 \ V 2 
r H = x - 1 = —-( — ) « 12,600 Mpc, (19) 

where the second equality assumes the standard values 
^m.o = 0.3 and fi V)0 = 0.7. This horizon distance r H is 
not the same as the particle horizon, but rather is essen- 
tially the Hubble radius. This distance scale rn provides 
an effective "boundary for microphysics" within the much 
larger space-time of the universe (for further discussion of 
horizons, see Kolb & Turner 1990 and Ellis & Rothman 
1993). 

For a flat universe with both matter and vacuum com- 
ponents, the scale factor a(t) increases according to the 
function 

'frd^)'"!-!^*]}" 1 

(20) 

At later times, the scale factor approaches the simpler 
asymptotic form 

a, 



>(*) = 



.,0 \l/3 



'4fi v ,o' 



explv 7 ^ H t] . (21) 



One can see immediately that the e-folding time scale of 
the future universe is r e = tt v ' H a w 17 Gyr. Fur- 
thermore, the scale factor a(t) approaches this asymptotic 
form on an even shorter time scale. If the full expression 
(eq. [20]) is written as the asymptotic form (eq. [21]) 
plus correction terms, those correction terms decay with 
a time scale r = (3-ffo-v/^v,o) _1 = T e /3 w 5.6 Gyr. This 
time scale — somewhat longer than the age of the solar 
system and appreciably younger than the current age of 
the universe — is a direct manifestation of the cosmologi- 
cal constant problem. 

Given that any extant cosmic structure has a sphere of 
influence (eq. [11]) and that the universe has a fixed hori- 
zon size Th, every structure will become isolated when the 
radius of its sphere of influence is stretched beyond the 
horizon, i.e., when a(t)r > Th- For a structure of mass 
Mobj , isolation occurs at a time t- lso given by 



4/ 3*(l 5 )3/2 



c 3 r fi 



GM obi O 



obj ^ ^m,0 



1/2 



(22) 



As a result, for vast majority of cosmological time, the 
cosmos will be divided into "island universes" in the sense 
that bound clusters of galaxies will retain their sizes (a few 
to several Mpc) while the distance between clusters grows 
exponentially (see also Chieuh & He 2002). Given typi- 
cal cluster sizes and separations, we predict that clusters 
will grow isolated in about 120 Gyr. For the particular 
values appropriate for the nearby Virgo cluster, equation 



(22) implies an isolation time of 132 Gyr (this time is the 
age of the universe at the time of isolation; since the uni- 
verse is already about 14 Gyr old, this event will occur 118 
Gyr from now). Structures with lower mass have smaller 
spheres of gravitational influence and require longer times 
to grow isolated. Our local group, with an estimated mass 
of M LG « (2.3 ± 0.6) x 10 12 M (van den Bergh 1999), 
will be isolated at cosmic age t = 175 Gyr. As another 
example, a star that is not gravitationally bound to a 
larger structure (e.g., one that has been scattered out of 
a galaxy) requires somewhat longer to become isolated — 
about 336 Gyr. For comparison, the lifetimes of the small- 
est, longest-lived stars are measured in tens of trillions of 
years (Laughlin, Bodcnheimer, & Adams 1997; hereafter 
LBA97), about one hundred times longer than the isola- 
tion time for galaxy clusters. For most of eternity, and 
indeed for most of the Stelliferous Era, clusters will be 
alone. Inside the galaxies, the expansion has essentially 
no effect, and star formation and stellar evolution con- 
tinue for trillions of years (AL97). When viewed on the 
large scale, however, these clusters will behave like point 
sources, pumping radiation into an ever-expanding void. 

4.2. Asymptotic structure of space-time 

In the long term, existing cosmic structures will remain 
bound, but will grow isolated. These structures will be 
embedded within an accelerating universe with a constant 
horizon scale (eq. 19). This process effectively divides the 
present-day universe into many smaller regions of space- 
time. These "island universes" display properties of a uni- 
versal nature. 

Every given "island universe" will approach a fixed over- 
density. The energy density contained with the horizon 
scale is equivalent to a mass Mh given by 

47T 



= ) r H = 
15 \3/2 



Since the mass contained within any given isolated clus- 
ter will be constant, the overdensity approaches a constant 
value. For our particular environment, the local group will 
remain bound with its mass of about Mlg = 2.3 xlO 12 
Mq. The mass in our local region is thus destined to be a 
minor perturbation on the cosmos itself, even within our 
local island universe. For cosmic ages older than 175 Gyr, 
the mass contribution to the universe contained within the 
local group is given by Soo — Mlq/Mh ~ 3 x 10~ 12 , only 
3 parts per trillion. 

Our numerical simulations indicate that cosmic struc- 
tures, from galaxies to clusters, tend to develop universal 
forms for their density profiles. As a result, every island 
universe will attain the same general form for its space- 
time. In particular, since the density profile attains the 
universal form described by equation (16), the line ele- 
ment ds 2 for the space-time within the horizon distance 
rn also attains a universal form. If we take the center of 
the coordinate system to be the center of the cluster and 
assume that the mass distribution is spherically symmet- 
ric, the line element can be written in the form 

ds 2 = -(l-A{r)- X 2 r 2 )dt 2 



+ {l - B(r) - xV) 1 dr 2 +r 



dfl 2 



(24) 
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where A(r) and B(r) depend on the mass distribution (see, 
e.g., Misner, Thorne, & Wheeler 1973). This form for the 
line element is that of a mass distribution embedded in de- 
Sitter space (see also Bardeen 1981, Mallet 1985, Chiueh 
& He 2002). As a result, there exists an outer horizon at 
r = x~ x ■ We provide a more detailed specification of the 
metric in a future work (Adams et al 2003). 

The outer horizon supports the emission of radiation 
through a Hawking-like mechanism (e.g., Fulling 1977, Bir- 
rcll & Davies 1982). As a result, the universe will be filled 
with a nearly thermal bath of radiation with characteristic 
wavelength A <~ rjj ~ x 1 ~ 12, 600 Mpc and character- 
istic temperature T ~ \ ~ A 2 /M pl ~ 10~ 33 cV ~ 10~ 29 
K. This bath of radiation will become the dominant back- 
ground radiation field at very late times (after about one 
trillion years - see the following section). 

4.3. Background radiation fields in an accelerating 
universe 

As the universe expands, all radiation fields are red- 
shifted to longer wavelengths. An important milestone 
is reached when the typical wavelength of a given radia- 
tion field grows longer than the cosmological horizon scale 
defined by equation (19). After this crossing, the pho- 
tons are larger than the largest "box" that the universe 
has to contain them. For later times, it no longer makes 
sense to describe the photons in terms of a distribution 
function. Inside the horizon, in the limit A ^> rjy, the 
background photons will appear as "slowly" varying elec- 
tric fields rather than as particles of light. The dominant 
background radiation field will be that produced by the 
horizon itself through a Hawking-like mechanism (see §4.2, 
Fulling 1977, Birrell & Davies 1982). 

Given the scale factor of the universe and the present 
day wavelength of a radiation field, it is straightforward to 
find the time at which the photons are stretched beyond 
the horizon scale, i.e., when Xa(t) > tr- For the cosmic 
background radiation, the present day wavelength (at the 
peak of the distribution) is about A = 0.1 cm and the 
photons cross the horizon at a time of 1120 Gyr. 

The cosmic background photons are stretched beyond 
the horizon well before the stars stop shining. Star forma- 
tion and stellar evolution will continue until the universe 
is tens of trillions of years old (AL97, LBA97). Suppose 
that stars continue to shine for 10 trillion years. By this 
late epoch, most of the remaining stars will be red dwarfs 
that emit light with a characteristic wavelength of A = 1 
fjbia = 10~ 4 cm. If this red light is emitted up to a time 
i* (« 10 12 yr) and observed at a later time t, its observed 
wavelength is given by 

A bs = Aemit ° bs = Aomit exp [(t - t*)/r e ] , (25) 

&emit 

where r e = 17 Gyr is the e- folding time of the future uni- 
verse. Starlight leaves the horizon when A bs > tr- Using 
this criterion in conjunction with equations (19) and (20), 
wc find that starlight is rcdshifted out of the horizon over a 
time interval of only At = (i — i*) = 1260 Gyr. Stellar evo- 
lution times — for the smallest stars — are much longer 
than the cosmological expansion times, so that photons 
are rapidly stretched beyond the horizon. Specifically, this 
stretching time is a small fraction of the Stellifcrous Era, 
the time over which the universe will contain substantial 
numbers of hydrogen burning stars, i.e., Ai/i* s=s 10~ 2 . 



4.4. Particle annihilation in an accelerating universe 

For material between galaxies — particles that are not 
bound to large structures — future evolution can continue 
through particle annihilation. The number density n of a 
given particle species is given by the evolution equation 



dn 

— +3Hn 
dt 



— (av)n 2 



(26) 



where (av) is the appropriate average of the interaction 
cross section and the relative velocity (see, e.g., Kolb & 
Turner 1990; see also Cirkovic & Samurovic 2001). For an 
accelerating universe with a cosmological constant, the so- 
lution to equation (26) can be found and written in terms 
of the scale factor, i.e., 



n(a) = n a a 3 < 1 



2r 



m,0 



n 



v,0 



On 



.00- 



-3 



(27) 

where no is the particle density at the present epoch and 
where we have defined T = (<rv)no/H . The leading fac- 
tor (a~ 3 ) represents the dilution of the number density due 
to cosmic expansion, whereas the second factor in brack- 
ets incorporates the effects of continued particle annihi- 
lation. In an accelerating universe, annihilation is highly 
suppressed. For example, in the asymptotic limit t — > oo, 
a — > oo , this factor becomes T = 1 + 2r( 1 — ^/f2 v .o) /3^ m ,o 
« 1 + 0.0048 [(crv)/barn-c]. For electron-positron annihi- 
lation, for example, the maximum correction term is less 
than a percent. For annihilation of cold dark matter par- 
ticles (thought to be the dominant matter contribution), 
the interaction cross sections typically lie in the range a ~ 
10~ 12 - 10~ 14 barn (e.g., Kolb & Turner 1990) and the 
speed v/c ~ 10 -3 . As a result, the already small correc- 
tion term (0.0048) is suppressed by an additional 16 orders 
of magnitude. This enormous suppression is driven by the 
relentless expansion of an accelerating universe. 

When the number density grows so diffuse that the uni- 
verse contains less than one particle per horizon volume, 
then individual particles are effectively isolated. The con- 
dition for such isolation can be written in the form 



, , /47T \V3 



(28) 



where the scale factor a is given by equation (20) and the 
horizon scale is given by equation (19). Adopting a present 
day number density of n = 10~ 6 cm~ 3 results in a particle 
isolation time scale of about 1060 Gyr. 

5. DISCUSSION AND CONCLUSIONS 

This paper explores the future evolution of a universe 
dominated by dark vacuum energy. Analytic estimates are 
compared to results of numerical simulations that follow 
the evolution of future structures in such an accelerating 
universe. 

For a universe with cosmological constant fi v ,o = 0.7, 
only those regions with present-day overdensities <5o > 17.6 
will remain gravitationally bound, in agreement with ear- 
lier estimates (Lokas & Hoffman 2002; Nagamine & Loeb 
2003). We generalize this result to include quintessence 
models with constant forms for the equation of state (Ap- 
pendix A). We have also derived the condition required for 
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test bodies to remain bound to existing structures (see eq. 
[11]) and verified its validity with numerical simulations 
(to within ~ 10%; see Figure 2). Any collapsed object - 
from a star to galaxy cluster — has a finite sphere of grav- 
itational influence in an accelerating universe, with radius 
r « 1 Mpc (Mobj/lO^Mo) 1 / 3 . For quintessence mod- 
els, we have derived an analogous result for the sphere of 
gravitational influence (see Appendix A and eq. [A7]). 

From a co-moving perspective, the large-scale appear- 
ance of the future universe is little changed from that of 
today. Matter in the cosmic web drains efficiently into 
collapsed halos that shrink in comoving coordinates. The 
halos are essentially frozen in place while their contrast rel- 
ative to the mean background grows with time. In physical 
coordinates, the view is rather different. The vast majority 
of the galaxies now visible are pulled out of the immediate 
horizon of any given bound structure (a cluster or group) . 
In the long term, only the cluster or group itself remains 
within the effective horizon scale of th = 12,600 Mpc. 

The long-term structure of space-time consists of a flat 
metric dimpled with isolated clusters that approach a fixed 
mass profile. We find that halo density profiles approach a 
form similar to, but steeper at large radii than, the NFW 
profile (equation [16]). It is important to emphasize that 
every halo grows isolated in the long term, i.e., every grav- 
itationally bound mass concentration ultimately becomes 
the only structure within its own island universe. In each 
such local region, the halo density takes the form shown 
in Figure 4 and the line clement of the space-time metric 
takes the form given by equation (24) . Although the halo 
mass varies from region to region, the form of the metric - 
and hence the geometry of space-time - is nearly univer- 
sal. In all cases, the halo mass in any region provides only 
a minor contribution to the overall mass/energy budget, 
with M ohj /M H ~ 10" 11 . 



As the universe continues to expand, and accelerate, 
cosmic radiation fields are redshifted to increasingly long 
wavelengths. After about one trillion years, the cos- 
mic background radiation (leftover from the big bang) is 
stretched beyond the horizon and the dominant radiation 
background is that emitted by the horizon itself through 
a Hawking-like mechanism. Many of the results of this 
investigation can be summarized in terms of the relevant 
time scales, which are listed in Table 1. To emphasize the 
mismatch between the various time scales, the table also 
lists the scale factor for each relevant epoch. For example, 
individual stars grow isolated in 336 Gyr (a = 2 x 10 8 ), 
but the longest-lived stars burn hydrogen for 17,000 Gyr 
(when a = 10 434 ). 

With the analysis complete, many of the time scales 
and length scales of the future universe can be understood 
in simpler terms. A dimensional analysis - presented in 
Appendix B - shows that the most important time scale 
is given by the asymptotic form for the Hubble parameter 
Hoo = y/QvfiHo ~ 59 km s _1 Mpc -1 . The accelerated 
expansion itself completely dominates the evolution of the 
universe as a whole, so that all of the time scales are de- 
termined by H^ 1 « 17 Gyr and logarithmic multiplying 
factors (see Appendix B). By comparison, time scales for 
stellar evolution (10 13 - 10 14 yr; AL97, LBA97) and dy- 
namical relaxation of galaxies (10 20 yr; BT87) arc much 
longer. 

The authors would like to thank Greg Laughlin for use- 
ful discussions. This work was supported by the Univer- 
sity of Michigan through a Regents Fellowship (MTB), 
a Physics Department Fellowship (RHW), by the Michi- 
gan Center for Theoretical Physics, and by NASA Astro- 
physics Theory Grant NAG5-8458. 



APPENDIX 

A. QUINTESSENCE 

In this Appendix we generalize our results to the case of a vacuum energy that depends on time, or equivalcntly, the 
scale factor a. For the sake of defimteness, we adopt the standard form for the vacuum energy equation of state, i.e., the 
vacuum pressure is given by 

Pvac = WPvac , (Al) 

where the parameter w is constant and lies in the range — 1 < w < 0. Current observations seem to indicate a somewhat 
smaller range — 1 < w ^ — 0.5 (e.g., Limin et al. 2000, Balbi et al. 2001). In fact, Spergel et al (2003) place a 95% 
confidence limit of w < —0.5 using a combination of the WMAP CMB data and the HST key project (Freedman et al. 
2001) value for the Hubble constant and find that w < —0.78 when additional constraints are added (from the SNIa- 
derived redshift distance relation, the 2dFGRS large-scale structure, and Lyman-a data, assuming a flat universe with 
constant w). For completeness, we consider here the full range of constant w values from to —1. With this equation of 
state, the scale factor evolves according to 

(l) 2 = H o{^m, a- 3 + n v . a-p} (A2) 

where the index p = 3(1 + w). 

The energy equation (6) that determines whether or not overdense regions collapse, and the fate of test bodies, can be 
written in the form 

(§) 2 = n v,o^ p ( 2 -0+ (fi m ,o + , (A3) 

where (3 measures the gravitational influence of an existing structure according to equation (7). The case of overdense 
regions can be considered by replacing (3 with f2 mi0 c>o, where So is the overdensity. 
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Table 1 

Time Scales and Scale Factors 



Event 


Time r (Gyr) 


aM 


Time scale for scale factor to approach exponential form 


5.6 




Inverse l~Tiil~)l~>le constant J~f ^ 


14 




e- folding time of the future universe {Hq^JVI^ o) _1 


17 




Current age of the universe 


13.7 


1 


Virgo Cluster leaves our horizon 


132 


1000 


The Local Group grows isolated 


180 


2 x 10 4 


Exiled stars become isolated 


336 


2 x 10 s 


Individual particles grow isolated 


1060 


6 x 10 26 


CBR photons stretch beyond the horizon 


1120 


2 x 10 28 


Optical photons stretch beyond the horizon 


1260 


10 32 


Lifetime of longest-lived stars 


17,000 


10 434 


End of the Stcllifcrous Era 


100,000 


^q2554 



In order to determine whether trajectories turn around (and hence remain bound or collapse), the right hand side of 
equation (A3) must vanish as before. In this case, however, the resulting cubic equation has time dependent coefficients 
and so the evolution of the scale factor must be considered simultaneously. As a result, we combine the two equations by 
changing the independent variable to a and thereby obtain 

W a 2 [n m , a- 3 + n vfi a~p] ' 1 ' 

Trajectories turn around when the right hand side of this equation vanishes. The minimum value of (3 required for such 
turnaround occurs when the right hand side of the equation has a double zero (both the right hand side and its derivative 
with respect to £ vanish). To find the critical value of (3, denoted here as /?*, we numerically integrate equation (A4) and 
iterate to find the value that provides a double zero. This procedure must be carried out for every value of w (or p). The 
resulting values of (3 are given below. The overdensities required for the collapse of future structures, for a cosmology 
with a given value of w, are given by 8q — /3*/^m,o- This quantity is shown in Figure A5. 
We also provide a simple fit to the numerical result: 

(3*(w) = /ft [1 + a(l + w) + 6(1 + wf + c(l + wf] , (A5) 

where /3g is the value for a cosmological constant ((3q « 5.3; §3.2) and where the coefficients are given by 

a = -2.33, 6 = 2.58, c=-1.20. (A6) 

This simple cubic fit reproduces the numerical results with an absolute error bounded by 0.035/3q ~ 0.19 (a relative error 
of a few percent - see Figure A5). A better fit could be obtained by using polynomials of higher order, but this level 
of accuracy should be adequate for most applications. With this fitting polynomial, the sphere of influence for existing 
structures, tq = (2G M ^j / (3* Hq ) 1 / 3 , can be written in the form 

r G « 0.7 Mpc (-Jg-) V Vo 2/3 [1 + «(1 + v) + 6(1 + ™? + c(l + ~ 1/3 • (A7) 

For accelerating universes, we can find the maximum distance that a light signal can propagate between now (the 
present epoch) and temporal infinity. This maximum distance r max is given by 



to 



cdt c f°° da 
o(t) = 



,2 



-1/2 

(A8) 



By solving the integral numerically and fitting the result, we can write the distance scale in the form 

r max H Q /c = I(w) w J [1 + 5(1 + w) + 6(1 + wf + 2(1 + w) 3 ] , (A9) 

where Iq = 1.141, d = 3.073, 6 = -12.39, and c = 37.13. This fitting function is valid for the range of equations of state 
— 1 < w < —1/2. As w — > —1/3, the integral (and hence r max ) becomes divergent. The result is shown in Figure A6. 

With the sphere of gravitational influence defined by equation (A7) and the maximum distance defined by equation 
(A9), we can define the isolation time ti so for structures through the relation a(t)rc > ^max- The asymptotic form for the 
scale factor a(t) is given by 

a(t)~ (|v^^ot) 2/P , (A10) 
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Fig. A5. — This plot shows the overdcnsity required for the 
collapse of future structures as a function of the parameter w ap- 
pearing in the equation of state for quintessence models. The 
sphere of gravitational influence for existing cosmological struc- 
tures is given by tq = (2GM ^,j / (3* Hq ) 1 ^ 3 , where the parameter 
f3* is related to the overdensity required for collapse of existing 
regions via f3* = 5o^lm,o (see text). The solid curve shows the nu- 
merically determined values; the dashed curve shows a cubic fit to 
the function (using eqs. [A5 — A7]). As w becomes more negative, 
the overdensity required for collapse becomes larger and the sphere 
of gravitational influence (for existing structures) grows smaller — 
the formation of future structure is more suppressed for smaller w. 




w 



Fig. A6. — The maximum distance that a light signal can propa- 
gate between the present epoch and temporal infinity for an accel- 
erating universe described by equation of state parameter w. The 
dashed curve shows a cubic fit to the numerically obtained result 
(see Appendix A and eq. [A9]). 



where p = 3(1 + w) as before. The isolation time is then given by 

i iso « -¥kLll( w )\j3*( w )]W c[2GM ohi H Q ]- 1/3 Y /2 . (All) 

Notice that this form applies only for values of p strictly greater than zero (w > — 1). To properly take the limit w — ► — 1, 
p — > 0, the function must include additional terms that are neglected in this approximation. 

B. DIMENSIONAL ANALYSIS 

In this Appendix, we present a dimensional analysis that illustrates the fundamental results of this paper in simpler 
terms. With the benefit of hindsight, we can conceptually reproduce many of the results of this paper. 

Because the universe is already dominated by its dark vacuum contribution, the future behavior of the universe is 
essentially one of exponential expansion at a well defined rate. This rate is set by the Hubble constant. For the sake of 
definitcness, we will use the asymptotic value of the Hubble constant = H ^^l v ,o ~ 59 km s _1 Mpc -1 . This rate 
also defines the basic time scale for the problem, i.e., r = H^ 1 = 17 Gyr. 

Using square brackets to denote the units of a given quantity (in terms of length L, time T, and mass M), we can list 
the variables that describe the the asymptotic universe via 

[ffoo] = T- 1 [G] = L 3 /MT 2 [c] = L/T . (Bl) 

These variables can be combined to produce a dimensionless field n if only if a mass scale M is introduced. The field 
n is then given by 

n .^. ,B2, 

If no additional variables are introduced into the problem - no additional entities are introduced into the universe - then 
typically Ho ~ 1, which in turn defines a mass scale for the universe. Inserting numerical values, we find Mq ~ 10 23 M Q 
(essentially the same result as that of equation [23]). 

Both our physical intuition and the results of our numerical simulations indicate that cosmic structure becomes frozen 
and bound astronomical objects - galaxies and clusters - grow isolated in the long term. The presence of a cluster or 
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galaxy introduces another variable into the problem, namely the mass scale M D bj. This scale, in turn, defines another 
dimensionless field II i given by 

H^^^^O-n, (B3) 

an incredibly small number. This quantity represents the overdensity of an isolated cluster embedded in an island universe 
at late times. 

Next, we want to define a length scale r associated with the galaxy or cluster itself. Such a length scale can be defined 
in several ways. The natural length scale of the universe is given by c/H^ = tr- Since the ratio tq/th is dimensionless, 
we expect that 

r = r H ^ = ^(^^Y, (B4) 

where the power-law index n is to be determined. If we argue that the length scale associated with the galaxy or cluster 
should be non-relativistic, then n must be chosen so that the scale ro does not depend on the speed of light. This constraint 
specifies n — 1/3 and hence 

r = (GMobj/ffo 2 ) 1 / 3 , (B5) 

which is the same as the gravitational sphere of influence defined by equation (11) (up to dimensionless factors of order 
unity) . Equation (B4) allows for a second "natural" length scale - that determined by eliminating the Hubble parameter 
by using n = 1. This choice results in the scale r = GM bj/c 2 , which is the length scale that determines the form of the 
functions A(r) and B(r) appearing in the metric (eq. [24]). 

Another result of this investigation is the time scales for which objects become isolated and radiation is stretched 
"beyond the horizon" . At this level of analysis, all of these time scales are the same and are determined by the asymptotic 
c-folding time t — H^ 1 w 17 Gyr. The more detailed mathematical analysis of the paper includes logarithmic correction 
factors, as listed in Table 1. For cosmological events, however, even the longest time scale is only 1260 Gyr or 74 e-folding 
times (it also turns out that ln[r^f/A] ~ 75 - for 'typical' astrophysical photons of wavelength A ~ lpm). For a universe 
with a cosmological constant, the basic result is that all future cosmological events must unfold with "nearly" the same 
time scale, given by H^ 1 . For comparison, stellar evolution time scales are determined by more complicated physics and 
span a wider range of time scales, both much shorter (~ 10 5 yr for star formation events; Adams & Fatuzzo 1996) and 
much longer (<~ 10 13 — 10 14 yr for the duration of the longest-lived stars; AL97, LBA97). Galactic evolution dynamical 
relaxation and evaporation - takes place over still longer times (~ 10 20 yr; Dyson 1979, Binney & Tremaine 1987, AL97). 
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